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Full jet reconstruction in heavy-ion collisions is a promising tool for the quantitative 
study of properties of the dense medium produced at RHIC. In addition to baseline 
measurements in p+p, results from d+Au collisions are needed to disentangle initial 
state nuclear effects from medium-induced fcx broadening and jet quenching. We present 
mid-rapidity inclusive jet px spectra and di-jet correlations (fcx) in 200 GeV p+p and 
d+Au collisions from the 2007-2008 RHIC run. 



1. Introduction 

Jets are coUimated sprays of particles produced in high energy particle collisions. 
They are remnants of hard-scattered partons, which are the fundamental objects of 
pQCD. At RHIC, jets can be used as a probe of the hot and dense matter created 

..19 

in heavy ion collisions .To quantify the signals observed in heavy ion collisions 
in comparison to p+p collisions, it is necessary to measure the cold nuclear matter 
effects in systems such as d-|-Au. 

Cold nuclear matter effects can be described by partonic rescattering ^ and 
by modification of parton distribution functions Si. Measurements of single particle 
RdAu have shown modification even at midrapidity^^. It is important to verify these 
observations using jets: they provide much higher kinematic reach as they are not 
prone to fragmentation biases and one can also avoid uncertainties coming from 
imprecise knowledge of fragmentation functions. An alternative approach to fulliet 
reconstruction is a method using di-hadron correlations, discussed in detail in^^. 

2. Jet reconstruction 

This analysis is based on ^snn = 200 GeV data from the STAR experiment, 
recorded during RHIC run 8 (2007-2008). The Beam Beam Counter detector, lo- 
cated in the Au nucleus fragmentation region, was used to select the 20% high- 
est multiplicity events in d-|-Au collisions. The Barrel Electromagnetic Calorimeter 
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(BEMC) detector was used to measure the neutral component of jets, and the Time 
Projection Chamber (TPC) detector was used to measure the charged component 
of jets. In the case of a TPC track pointing to a BEMC tower, its momentum was 
subtracted from the tower energy to avoid double counting (electrons, MIP and 
possible hadronic showers in the BEMC). 

To minimize the effect of BEMC backgrounds and dead areas, the jet neutral 
energy fraction is required to be within (0.1,0.9). An upper < 15 GeV/c cut 
was applied to TPC tracks due to uncertainties in TPC tracking performance at 
high-pT in run 8 (under further investigation). The acceptance of TPC and BEMC 
together with experimental details (calibration, primary vertex position cuts) limit 
the jet fiducial acceptance to \r]\ < 0.55 [R — 0.4), \ri\ < 0.4 {R = 0.5), where R is 
the resolution parameter used in jet finding. 

Recombination jet algorithms kt and anti-kt, part of the FastJet package 
are used for jet reconstruction. To subtract d+Au underlying event background, a 
method based on active jet areas ^ is applied event-by-event: p^:'^'^ = pCandidate _ 
p ■ A, with p estimating the background density per event and A being the jet 
active area. Due to the asymmetry of the colliding d-f Au system, the background 
is asymmetric in 77. This dependence was fit with a linear function in rj and included 
in the background subtraction procedure. 

Pythia 6.410 and GEANT detector simulations were used to correct for exper- 
imental effects. Jet reconstruction was run at MC hadron level (PyMC) and at 
detector level (PyGe). To study residual effects of the d-l-Au background (such as 
background fiuctuations), a sample with added background (PyBg) was created by 
mixing Pythia events with 0-20% highest multiplicity d-|-Au events (minimum bias 
online trigger). This mixing was done at detector level (reconstructed TPC tracks 
and BEMC towers). 



3. Nuclear broadening 

Comparing azimuthal correlations of jets in di-jet events in p-l-p and d+Au can 
provide information on nuclear fcx broadening. To increase di-jet yield, BEMC high 
tower (HT) online trigger was employed (one tower with Et > 4.3 GeV) for both 
p-|-p and d+Au data (run 8). Resolution parameter R = 0.5 was used for jet finding 
and a cut px > 0.5 GeV /c applied for tracks and towers to reduce background. 

To select a clean di-jet sample two highest energy jets (pt,i > Pt,2) in each 
event were used, with pT.2 > 10 GeV/cl^. Distributions of fcx^raw = PT,isin(A<^) 
were constructed for di-jets and Gaussian widths, (Xfe^ , were obtained for the two 
jet algorithms and two (10 - 20 GeV/c, 20 - 30 GeV/c) pT,2 bins. 

Detector and residual background effects on the fcx widths were studied by 
comparing PyMC, PyGe and PyBg distributions as shown in Fig. [TJ The widths 
are the same within statistical errors, most likely due to the interplay between jet 
Pt and A0 resolutions. Fig.[2]shows an example of the fcx.raw distributions for data. 
The Gaussian fit to p-f p data is not ideal and the precise shape of the distribution 
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is under study. RMS widths of these distributions have therefore been checked and 
they agree with the sigma widths of the fits. The values extracted from the Gaussian 
fits are crf+P^^ = 2.8±0.1 (stat) GeV/c and cr^+f^ = 3.0±0.1 (stat) GeV/c. Possible 
nuclear fcx broadening therefore seems rather small. 

The systematic uncertainties on extracted fcx widths come from Jet Energy 
Scale (JES) uncertainty (discussed in more detail in section [5]) and from the way 
they were extracted. The latter is estimated to be 0.2 GeV/c and includes a weak 
dependence on the |A(/)— 7r| cut for back-to-back di-jet selection (varied between 0.5 
and 1.0), differences in thepx,2 range and jet algorithm selections and the precision 
with which the detector and background effects were found to be negligible. 



PyMC STAR Preliminary 

\Ss = 200 GeV, anti-kt, R = 0.5 
o„ = 2.1 1± 0.11 GeV/c 
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^*sin(A(|))[GeV/c] 

Fig. 1. Distributions of feT.ra 



PyGe STAR Preliminary 
Vs = 200 GeV, anti-kt, R = 0.5 
<'.t ™ = 2.19+ 0.10 GeV/c 
.+ 




PyBg STAR Preliminary 

s/s = 200 GeV, anti-kt, R = 0.5 
a„„, =2.1810.09 GeV/c 



p^^*sin(A ^) [GeV/c] 
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: Pt,i sin(A0) for simulation (10 < Pt,2 < 20 GcV/c). 



4. Inclusive jet spectra 

Run 8 d-f Au data with a minimum bias online trigger were used for this study. lOM 
0-20% highest multiplicity events after event cuts were used for jet finding (anti-kt 
algorithm) with resolution parameter R = 0.4, pT > 0.2 GeV/c cut was applied 
to tracks and towers. The jet pt spectrum is normalized per event and the high 
multiplicity of d-|-Au events also guarantees the trigger efficiency is independent of 
the Pt of the hard scattering. Therefore, no correction related to trigger is applied 
to jet Pt spectra. 

A bin-by-bin correction is used to correct jet pt spectrum to hadron level. It 
is based on the generalized efficiency, constructed as the ratio of PyMC to PyBg 
jet Pt spectra, applied to the measured jet pt spectrum. It therefore corrects for 
detector effects (tracking efficiency, unobserved neutral energy, jet pT resolution) 
as well as for residual background effects. As the impact of these effects on jet pt 
spectrum differ substantially depending on jet pt spectrum shape, the shapes have 
to be consistent between the PyBg and the measured jet pT spectra. Fig. [3] shows 
that this is indeed the case. 

As the jet pT spectrum is very sensitive to the jet energy scale, an additional 
correction was applied here to account for the lower TPC tracking efficiency in 
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d+Au compared to that from the used p+p Pythia simulation. The d+Au efficiency 
was determined by simulating single pions and embedding them at the raw detector 
level into real d+Au minimum bias events. The tracking efhciency in the Pythia 
simulation was then artificially lowered, prior to jet finding at PyGe and PyBg 
level, so that it matches the one obtained from d+Au embedding. 

To compare per event jet yield in d+Au to jet cross section in p+p collisions, 
an input from MC Glauber study is utilized: (iVbin) = 14.6 ± 1.7 for 0-20% highest 
multiplicity d+Au collisions and crinei.pp = 42 mb. These factors were used to 
scale the p+p jet cross section measured previously by the STAR collaboration H^l 
using Mid Point Cone (MPC) jet algorithm with R = 0.4. The resulting d+Au 
jet pt spectrum is shown in Fig. 0] together with the scaled p+p jet spectrum. 
The systematic errors are indicated by dashed lines and by the gray boxes. The 
dominant contribution to p+p systematic uncertainty is the Jet Energy Scale (JES) 
uncertainty. Within these systematic uncertainties, the d+Au jet spectrum shows 
no significant deviation from the scaled p+p spectrum. 

5. Systematic uncertainties 

The JES uncertainty dominates the uncertainties of d+Au measurement and is 
marked by the dashed lines in Fig. ID Part of it comes from the BEMC calibration 
uncertainty of 5%: it's applied to the neutral component of the jet (« 40% for 
minimum-bias online trigger, w 55% for HT online trigger). 

An uncertainty of 10% in TPC tracking efficiency is applied to the charged 
component of jets. Embedding of jets into real d+Au events at raw detector level 
will allow to decrease this uncertainty in the future. As JES uncertainty is expected 
to be largely correlated between run 8 p+p and d+Au data, we plan to measure 
jet pt spectrum in run 8 p+p collisions to decrease uncertainties in i?dAu- 

Caution is needed due to the use of different pseudorapidity acceptances and 
different jet algorithms in Fig. |4l The effect of ry acceptance on jet pt spectrum is 
illustrated in Fig. [5] using Pythia simulation: it is less than 10% in the pt range 
covered by the present measurement. The effect of jet al gori thm is mainly due to 
hadronization and can exceed 20% for pT < 20 GeV/cM Therefore, the same 
acceptance and the same algorithm have to be used in p+p and d+Au to obtain 

jet i?dAu- 

6. Summary 

Di-jet fcx widths were measured in 200 GeV p+p and d+Au collisions: cr^^^ = 
2.8 + 0.1 (stat) GeV/c, cr^Jf^" = 3.0 + 0.1 (stat) GeV/c. No significant broadening 
due to Cold Nuclear Matter effects was observed. Jet px spectrum from minimum 
bias 200 GeV d+Au collisions is consistent with the scaled p+p jet spectrum within 
systematic uncertainties. Precise tracking efficiency determination from jet embed- 
ding in raw d+Au data and jet cross section measurement in run 8 p+p data will 
allow to construct jet i?dAu- 
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p+p STAR Preliminary 

. \/s = 200 GeV, anti-kt, R = 0.5 
= 2.83 ±0.05 GeV/c 
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' sin(A (])) [GeV/c] 



d+Au STAR Preliminary 

\[s^ - 200 GeV, anti-kt, R = 0.5 
2.98 ± 0.08 GeV/c 




-2 2 4 6 8 

sin(A (|)) [GeV/c] 



^T,1 > T/ 1. J 1-.^^ 

Fig. 2. Distributions of fcx,raw for p+p, d+Au (10 < Pt,2 < 
20 GeV/c). 
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Fig. 3. Ratio of jet pT spectra 
between d+Au and simulation. 
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Fig. 4. Jet 

p+pHni. 
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Pt spectrum; d+Au collisions compared to scaled Fig. 5. Effect of r} acceptance on jet 

Pt spectra (Pythia). 
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